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The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic 
energy are computed using direct simulation data of a turbulent channel flow. The budget 
data reveal that all the terms in the budget become important close to the wall. For 
inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return 
term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The 
rapid and return terms play different roles depending on the component of the term. A 
split of the velocity pressure-gradient term into a redistributive term and a diffusion term 
is proposed, which should be simpler to model. The budget data is used to test existing 
closure models for the pressure-strain term, the dissipation rate, and the transport rate. 
In general, further work is needed to improve the models. 

1. Introduction 

The advancement of large-scale computers has led to the wide use of turbulence 
models to predict turbulent flows. At the 1980-81 AFOSR-HTTM-Stanford conference, 
phenomenological models of all categories (one-, two-, and Reynolds-stress-equations) 
were used to model various flows. The anticipated result, that transport models for the 
Reynolds-stress equations will outperform eddy-viscosity-type closures, did not happen. 
In fact, in simple flows, mixing-length models produced results as good or better than 
Reynolds-stress models. However, the same Reynolds-stress models were used in a variety 
of flows to simulate quantities that simpler models could not model. The possibility of 
using the same model with the same constants for a variety of flows is the main motivation 
behind developing models at the Reynolds-stress level. 

According to Lumley (1978), the history of “higher-order” modeling dates as far back 
as Kolmogorov (1942), who was the first to suggest the characterization of turbulence by 
its intensity and scale. Chou (1945 a) considered the full Reynolds-stress equations and 
the triple-correlation equations to close the averaged Navier-Stokes equation. He suggested 
that the use of equations up to third moments is sufficient to characterize turbulent flows 
and, with simplifications, he predicted the mean velocity profile of a turbulent channel 
flow (Chou, 1945 b) which is the flow of interest in this work. 

Rotta (1951a) developed closure models for the Reynolds-stress equations and ad- 
vanced models for the pressure-strain and the dissipation rate terms that are the basis 
for several of the present day models. In a followup paper, Rotta (1951 b) introduced an 
equation for the turbulence length scale and tested his models by simulating a turbulent 
channel flow. He found that the peak in the turbulence intensity could not be predicted 
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by his theory. Davydov (1959) proposed closures at the Reynolds-stress and the triple- 
correlation equations levels. Later, Davydov (1961) proposed closures for the equation of 
the dissipation rate of the turbulence kinetic energy (e) which are used by most present-day 
c-equations. Daly & Harlow (1970) (hereafter referred to as DH) used various ideas to close 
the Reynolds-stress equations and used gradient-transport to model the turbulent diffusion 
terms. They added the scalar equation representing the trace of the dissipation rate tensor 
to the Reynolds-stress equations. They also computed the channel flow with mixed success 
in predicting the Reynolds-stresses. Launder, Reece, & Rodi (1975) (hereafter referred to 
as LRR) advanced closure forms to the Reynolds-stress equations that were tested for a 
variety of turbulent flows. Their models were developed for high Reynolds number flows. 
Hanjalic & Launder (1976) (hereafter referred to as HL2) extended the model of LRR 
to simulate low-Reynolds-number turbulence and the near-wall region. HL2 computed 
the channel flow and found that their model underpredicts the peak in the turbulence 
kinetic energy by 30%. More recently, Shih & Lumley (1986) proposed closure forms to 
the pressure-strain terms that seem to hold promise but did not use the model to compute 
to the wall (wall functions were used in the near-wall region). 

All of the work on turbulence modeling development uses indirect methods to test 
the various closure models. Because of the difficulty in measuring pressure and velocity 
with sufficient spatial accuracy, direct comparison of the closure models with experimental 
data has not been possible. Often, the adequacy of the model is judged by computing a 
flow with the model and by comparing the predicted mean velocity and Reynolds stresses 
with experimental data. With the advent of large-scale computers and new algorithm 
developments, direct simulations of turbulent flows are now possible at moderate Reynolds 
numbers. These simulations are being used to compute the terms in the budget of the 
Reynolds stresses. It is then possible to test the closure models by direct comparison of 
the closure formula with the term being modeled. 

A comprehensive testing of one point closure models using simulation data was carried 
out by Rogallo (1981) where he tested different pressure-strain models for homogeneous 
flows under mean strain and shear. Recently, direct simulations of simple inhomogeneous 
flows have been carried out by Kim, Moin & Moser (1987) (hereafter referred to as KMM), 
Moser & Moin (1984), and Spalart (1986a, 1986b). KMM computed the fully turbulent 
channel flow using two different grid resolutions (2 x 10 6 and 4 x 10 6 grid points) and 
found no dependency of the results on grid resolution. They compared their results with 
experimental data and were able to study the behavior of turbulence correlations near 
the wall. Although there exist some disagreements between experimental data and com- 
puted results, especially in the near wall region, the overall agreements were good. KMM 
attributed the differences to possible probe errors due to wall proximity. Moser & Moin 
computed a channel flow with mild streamwise curvature. They showed that the presence 
of the Gortler vortices contribute substantially to the shear stress. They also computed the 
terms in the budget of the Reynolds stresses. Spalart (1986a) computed boundary-layer 
flows with favorable pressure gradient and more recently (Spalart, 1986b) computed tur- 
bulent boundary layers with zero pressure gradient up to Reg = 1410. These simulations 
were carried out at moderate Reynolds numbers and were used to compute the Reynolds 


2 



stresses and the terms in their budgets. 

The objective of this work is to present the budget data for the Reynolds stresses 
and the budget data for the dissipation rate of the turbulence kinetic energy using the 
flow fields of KMM. These budgets are needed for model development because they will 
provide a direct means to evaluate closure models. The budget data for the dissipation 
rate of the turbulence kinetic energy have eluded measurement techniques and will provide 
valuable guidelines for model developers and model testing. We will use these budgets to 
evaluate some existing algebraic models for the dissipation rate of the Reynolds stresses. 
Closure models for the budget of the dissipation rate of the turbulence kinetic energy will 
be evaluated. We will also compare with the data some existing models for the pressure- 
strain and the turbulent transport rate of the Reynolds stresses. 

2. Channel data analyses 

The use of direct simulation to compute turbulent flow fields will be restricted to 
simple flows for the forseeable future. The flows of engineering interest are simulated using 
the averaged Navier-Stokes equations in conjunction with closure models. With the advent 
of large computers, the trend will be towards using closures at the level of the budget for 
the Reynolds stresses. 

2.1. Reynolds stress budget 

The transport equations for the Reynolds stresses are derived from the Navier-Stokes 
equations by ensemble averaging the equations, then deriving equations for the fluctuat- 
ing stresses and ensemble averaging these equations. For incompressible turbulent flow, 
the transport equations nondimensionalized with u\/u (the wall-shear velocity, u T = 
y/vU, y | wall’ an ^ ^he kinematic viscosity, u) are given by 

~£^ u i u i = Pi) -h Pi) T Dij + II jj — (1) 

where D/Dt = d/dt + U k d/dx k , and the terms on the right-hand side of the above equation 
are identified as follows: 


Production rate 


Dissipation rate 
Turbulent transport rate 


Viscous diffusion rate 


Velocity pressure-gradient term 


Pii=~ 


u 'i u 'k U hk + u ') U 'k U i,k 


€ij — 2 u' ik u'- k 


Tij = -(«•«'<),* 
Dij =(u' i u' j ) <kk 
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Repeated indices imply summation over 1,2,3 and the indices (1,2,3) are used to denote 
the streamwise, x + , normal to the wall, y + , and spanwise, z + , directions respectively. In 
the above equation and in what follows, p' is a non-dimensional kinematic pressure. All 
quantities are { } + quantities, but the superscript + will be used for the coordinate vari- 
ables only to simplify the notation. In a fully developed channel, the flow is homogeneous 
in the streamwise ( x* ) a nd th e spa nwise (z + ) directions. The relevant non-zero stresses in 
this case are u\u\ y u' 2 u 2 , u' 3 u 3 , and u\u 2 . Figures 1-4 show the terms in the budget of these 
stresses using the flowfields of KMM. The simulation flowfields are for a channel flow at 
Reynolds number Re T = u T 6/u = 180 based on wall-shear velocity (tt T ) and channel half 
width (6). This corresponds to Reynolds number 3200 ba sed o n mean -c enter line velocity 
and 6. The budget for the turbulence kinetic energy k = |(u'jtt'j +u 2 u' 2 + u 3 u 3 ) is shown in 
Figure 5. The profiles for the different terms in these budgets (scaled by u\/u) are similar 
in shape and magnitude to those of Moser & Moin (1984), and Spalart (1986b). 

To analyze the near-wall asymptotic behavior of the different terms in the budget 
equations, we will expand the instantaneous velocity and pressure in Taylor series about- 
the-wall values as follows: 


«1 = &il/ + +ci2/ + +d,!y + + ... 

u ' 2 = c 2 y + +d 2 y + + ... , 9 . 

i , -l +2 , +3 

«3 = &3 y +C 3 y +d 3 y + ... 

p' = a p +b p y + +Cpt/ + +d p y + + ... 

where the coefficients « p , 6 j,& 2> ... are functions of « + , z + , and t. The coefficient c 2 in the 
u 2 expansion is related to the coefficients b\ and b 3 through the continuity equation 

2c 2 = —(f>i,i + 63 , 3 ) (3) 


The first coefficient, a p , in the pressure expansion is related to the coefficients cj and C3 
through the x- and z-momentum equation 


a p,i = 2 ci 

«p, 3 = 2 c 3 


which implies that 


Ci ,3 = c 3) j 


(5) 


The second coefficient, 6 p , in the pressure expansion is related to the coefficient C2 through 
the y-momentum equation 

6 p = 2c 2 (6) 

Finally, the third coefficients, d 3 and d 3 , in the velocity expansions are related to and 
63 , y-partial derivative of the x- and z-momentum equations 


6di — 6i,t — 2&i t n — 63,31 — 61,33 

6d 3 = 63, t — 263,33 — &i,i3 ~ ^3,11 
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Using Equations ( 2 ) and (3) in the expression for the Reynolds shear stress, —u\v 2 , yields 


—u\u' 2 


6j6 


3,3 


-v 


+ 


( 9 ) 


Chapman k Kuhn (1986) have used the “splatting argument” to show, and the data of 
KMM confirm, that near the wall the first term in the expansion of u \u' 2 is positively 
correlated. The data of KMM show that 6163,3 = 1.4 x 10~ 3 . In the fully developed 
channel, the mean velocity near the wall will vary as, 

U = y + — y +2 /(2fie r ) + ... (10) 


Substituting Equations (9) and (10) into the expression for the production rate term 
of the u\ u'j yields 

-Pi 1 = 6163,3 y + + ... (11) 

Figure 1 shows that the production term in the u\u\ budget is the dominant “producing” 
(positive) term in the range t/ + > 10. Moving towards the wall, the turbulent transport 
term becomes positive at around y + ss 10. Then the viscous diffusion term becomes 
positive at y + % 5 , while the production term *->• 0 as . Taylor-series expansion about 
the wall of the viscous diffusion rate term yields 

D 11 = 26361 T 126iCij/"*" 4- ... (12) 

where from the data 6161 = 0.13. At large y + the dissipation rate term and the velocity 
pressure-gradient term are both negative, and of the same order of magnitude, and they 
balance the production rate. Moving towards the wall the velocity pressure-gradient term 
tends towards zero. The asymptotic behavior close to the wall of II n is 

nn = -46^y + + ... (13) 

where from the simulation we have — 6jC] % 8.5 x 10 -4 . The dissipation rate term remains 
large at all y + and does not vanish at the wall 

en = 26 i 6 i -f 861 ci y + + ... ( 14 ) 

The above expansions show that close to the wall the dissipation rate balances the viscous 
diffusion rate plus the velocity pressure-gradient term. At the wall the viscous diffusion 
rate is equal to the dissipation rate. 

The budget for the u 2 u 2 component of the Reynolds-stress tensor (fig. 2 ) shows that 
the turbulent transport rate T 22 is of the same order as the other terms through most of 
the channel. Close to the wall this term decays faster than the other terms. The u' 2 u 2 
budget does not have a production term but the velocity pressure-gradient term is the 
dominant producing term. The dissipation rate term is the dominant consuming term. 
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The viscous diffusion rate term is small compared to the other terms except very close to 
the wall (y + < 20). Taylor-series expansion of the viscous diffusion rate term, £> 22 , in the 
near- wall region yields 

D 22 = 12c 2 c 2 y +2 -I- ... (15) 

From the simulation we have c^ci — 7.4 x 10“ 5 . Equation (15) shows that close to the 
wall the diffusion rate is positive. Expansion of the dissipation rate term e 22 yields 

e 22 = 8 c 2 c 2 y + -f ... (16) 

Finally, expansion of the velocity pressure-gradient term yields 

n 22 = -4c^y +2 + ... (17) 

As expected, close to the wall II 22 balances with Z) 22 — c 22 . The expansion, Equation (17), 
shows that II 22 is negative close to the wall. 

The budget for the u' 3 u ' 3 component of the Reynolds-stress tensor (fig. 3) shows that 
away from the wall the “producing” term is the velocity pressure-gradient term. The 
dominant consuming term is the dissipation rate term. Moving towards the wall, the 
velocity pressure-gradient term decreases as 


II33 = — 4&3C3j/ + + ... (18) 

From the simulation we have — 63C3 % 4.4 x 10~ 3 . The viscous diffusion rate term becomes 
important close to the wall and reaches a maximum at the wall 


D 33 — 26363 + 126 3 c 3 y + + ... (19) 

From the simulation data we have 6363 = 3.76 x 10 -2 . The dissipation rate term also 
reaches a maximum at the wall. The near-wall behavior of c 33 is given by 


C33 = 26363 + 86 3 c 3 y + + ... (20) 

At all y + the turbulent transport rate T 33 remains small compared to the other terms. 
This is in contrast with T 22 which is of the same order as the other terms in the budget of 
u 2 u 2 . 

As in the budget of u'jUj , the budget for u\u 2 (Fig. 4) is dominated by P 12 (production 
of — u\u 2 ). Away from the wall, the velocity pressure-gradient term balances with P J2 , 
while the other terms are small. Moving towards the wall, the turbulent transport term 
becomes important. Very close to the wall, the dissipation rate term and the viscous 
diffusion rate term become important. The sum of the two viscous terms (the viscous 
diffusion term and the dissipation rate term) yields a term that is small throughout the 
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channel, indicating that the viscosity plays a minor role in the dynamics of the Reynolds 
shear stress. The asymptotic behavior of the various terms close to the wall is given by 


P 12 = ~c^cjy ~ f4 + ... 
ti2 = 4&jc 2 y + + ••• 

T \2 = — 56ic 2 c 2 y + * + ... (21) 

D 12 =6fei c 2 y + + ... 

H 3 2 =— 26 ic 2 y" 1 ’ + ... 


where £>ic 2 % — 7. x 10 4 . The above expansions show that P 12 and Ti 2 decay much faster 
than the other terms as y + 0. 

The budget for k (fig. 5) is one half the sum of the budgets of the diagonal components 
of the Reynolds- stress tensor. We point out that away from the wall (y + > 30) the 
argument often used that P*. = e holds acceptably well. Moving towards the wall, the 
turbulent transport rate becomes important. It is a consuming term in the 30 > y + > 8 
range and a producing term very near the wall. In effect it is transporting energy from 
the maximum source towards the wall. As we get closer to the wall, the dissipation rate 
balances with the viscous diffusion rate plus the pressure diffusion rate. At the wall, the 
dissipation rate is non-zero and is equal to the viscous diffusion rate. The data show that 

(\waii = 0.166 

If we compare the above budget to the budget of Laufer (see, for example, Townsend, 1976), 
we find that both the turbulent transport rate and pressure diffusion rate were overesti- 
mated by the measurements. The viscous diffusion rate is underestimated by Laufer’s 
data, which yields a lower dissipation rate at the wall. It should be noted that there is a 
large difference between the Reynolds numbers of the simulation and Laufer’s data. We 
expect the turbulent transport terms which are large scale dependent terms, to be less 
sensitive to Reynolds number dependence than the dissipation rate terms, which are small 
scale dependent terms. Our data is consistent with the data of Moser & Moin (1984) for a 
flow in a curved channel, and Spalart (1986a & b) for flows over a flat, plate. Near the wall, 
all simulation data show that the pressure diffusion rate term remains small compared to 
the other terms. However, very close to the wall the pressure diffusion rate term is of the 
same order as the difference between the dissipation rate and the viscous diffusion rate. 

2.2. Dissipation rate budget 

In the discussion of the Reynolds- stress transport, we have seen that far from the 
wall the dominant terms are the production rate, dissipation rate, and velocity pressure- 
gradient terms. The production rate term is a function of the Reynolds stresses and 
the mean velocity, and it does not need modeling. However, the rest of the terms need 
modeling. 

A set of equations describing the dynamics of can be derived from the Navier- 
Stokes equation, but doing so will introduce six more equations and more terms to be 
modeled. The alternative used in the Reynolds-stress modeling is to model the dissipation 
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rate tensor in terms of the Reynolds stresses and a turbulence time scale. In the case of 
isotropic turbulence, the combination of the turbulence kinetic energy and the dissipation 
rate term provides a time scale for the decay of the turbulence. The approach suggested 
by Davydov (1961), and taken by DH and Hanjalic & Launder ( 1972) (hereafter HL1) is 
to introduce the scalar trace of the dissipation rate tensor 


e = 



( 22 ) 


and to model in terms of u^u'- and t. In this way one equation describing the evolution 
of e is needed for closure. The equation for e derived from the Navier-Stokes equation is 
given by _ 

JU = P * + P] + P\ + P? + T t + n« + D ( - T ( 23 ) 

We can identify the different terms on the right hand side as (rate of ...) 


p € = - 2u i,j u k,i S * 


P* yk u\ fm Skm 


P f ~ 2 U 'k U 'i,mUi,kT 

P? =-2u! ,u! 


l i,k u i,m k,m 


Tf = -(K U 'i,m U 'i,m),k 


=--(/>'m<, m ),fc 


D* =e 


y kk 




Production by mean velocity gradient 
Mixed production 
Gradient production 
Turbulent production 
Turbulent transport 

Pressure transport 
Viscous diffusion 
Dissipation 


where Sij = {Uij + Uj,i)/2 is the mean strain-rate. Tennekes & Lumley (1972) analyzed 
the vorticity fluctuation budget, which is related to the above budget for homogeneous 
turbulent flows. They inferred from an order-of-magnitude analysis that, in the high- 
Reynolds-number regime, the turbulent-production rate (P t 4 ) and the dissipation rate (T) 
dominate the balance equation. However, they point out that the difference of these terms 
yields a term of the same order as the other terms. The various terms in the balance 
equation for e are shown in figure 6. The errors in the budgets are expected to be highest 
for this case, because the computation of the terms in Eq. (23) involve correlations of the 
derivatives. There are two sources of errors, numerical (differencing errors) and statistical 
(limited sample size). Unfortunately, we cannot assess these errors directly, for example 
by running a finner grid for longer time. But an indirect measure of these errors is the 
imbalance in the budget. In the case of the dissipation budget, the imbalance in the budget 
is less than 2% (of T at the wall) throughout the channel (less than 0.04% for y + > 20). 
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The present results indicate that P f and T are the dominant terms in the core region 
of the channel, in agreement with Tennekes and Lumley’s analysis. Near the wall however, 
these terms remain large but are not larger than the other terms. Close to the wall (y + < 8) 
the production rate P' becomes of the same order as P *. In the range 6 < y + < 15 the 
mixed-production rate (P 2 ) is of the same order as P*. 

Very close to the wall (y + < 4 ), the profile of T goes through a local minimum at 
y+ 2 and reaches a maximum at the wall. This behavior is an artifact of the splitting 
of the u' u' .. term into a viscous diffusion rate and a dissipation rate term. Figure 7 

t f 771 X 1 771 AC AC 

shows that — T + D ( does not exhibit this peak and shows a monotonic decrease towards 
the wall. The total term might be simpler to model. The models for the different terms 
in the e-equation will be discussed in the next section. 

Taylor-series expansion of the terms in the e-equation yields the following: 


Pi =2fc 1 fe 3 ,32/ + + 0(j/ +3 ) 

P; - 4 77 ^y+ 2 + 0 (y +3 ) 

p l + °(y +3 ) 

Pf =3(616163,3 + 636361,1 )y + + 0(y + ) 

T ( = (616163,3 + 636361,1 )y + + 0 (y+ ) 

n f = — 8 c^ + Q(y- f ) 

D ( =4 {6] ,1 6] ,1 + 61,361,3 + 63,363,3 + 63,163,1 + 2 (ciC] + c 2 c 2 + c 3 c 3 ) + 2 c 2 c 2 } + 0 (y + ) 

T = 4 {61,161,] + 6], 361, 3 + 63,363,3 + 63,163,1 + 2 (c]Ci -|- c 2 c 2 + C3C3)} + 0 (y + ) 

( 24 ) 

At the wall the pressure transport rate term balances with the diffusion rate term and 
the dissipation rate term. P t 3 is of the order l/Pe T and can be neglected relative to the 
other terms. The above expansions show that P 2 = 0 (y +2 ), and the other production 
terms (P t J and P t 4 ) are 0 (y + ) close to the wall. 

2.3 Velocity pressure-gradient terms 

We have seen in the budget of the Reynolds stresses that the velocity pressure-gradient 
terms play the dominant role in energy redistribution among the components. The velocity 
pressure-gradient term expression has been extensively analyzed and modeled by various 
groups. Unfortunately, because of the lack of experimental data, the present models have 
not been tested. The usual approach in the analysis has been to split the expression into 
a redistributive part and a transport part. The redistributive part is then analyzed with 
the aid of the pressure Poisson equation. 

2 . 3 . 1 . Pressure-strain split 

The velocity pressure-gradient term is usually split into a pressure transport term and 
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a redistributive term 


Rij = ~ (p'{u[6j k + u'jSik)) k + <j>ij 


(25) 


where <fijj = 2p'sij (and = ( u\ ■ + u' ti )/2) is known as the pressure-strain term. Sub- 
stituting the Taylor series expansion for pressure and velocities in the definition of <f>ij 
yields 

<f>u = 2a p b iti y + + ... 


4 > 22 = 4 0pC 2 t/ + + ... 
4> 33 = 2 a p bj ^y + + 
4> 12 = + ... 


(26) 


For the case of a fully turbulent channel flow, the flow is homogeneous in the streamwise 
and spanwise directions and the split is irrelevant in these directions; IIu and II33 are the 
pressure-strain terms. The Taylor series expansions (eq. 26) show that the diagonal terms 
of <t>ij 1— > 0 as y + t— » 0, while the </> 12 component does not vanish at the wall. 

Figure 8 shows the distribution across the channel of the terms obtained by splitting 
II 22 into a transport term and a pressure-strain term. In this case the split produces 
terms that highly emphasize the presence of the wall. The reversal of the sign of the 
pressure-strain term near the wall was attributed by Moin & Kim (1982) to the “splatting” 
effect. The split chosen in equation (25) is not unique (Lumley, 1975); Launder (private 
communication, 1985) pointed out that the use of a decomposition suggested by Lumley: 


+ 4>ii 




(27) 


where 


•Kj = - ((p , (“i*j* + Uj<i*)) k - + p , Wj+*y 


(28) 


would reduce the negative levels of <f > near the wall to about one third their current 
levels. A split suggested by the balance equation for the anisotropy tensor of the Reynolds 
stresses is as follows: 


u\u\ 


where 


II ij — <f>ij “b 

2/> 

(29) 

1 

'f* 

II 

u i u 'i \ 

2k n “) 

(30) 


Note that the trace of is zero and therefore (^,j is redistributive. Figures 9-12 show the 
distribution across the channel of the terms in Eqs. (29) and (30). The negative levels of 
the redistributive 22 component near the wall are substantially reduced. In addition, the 
split suggests that a model for the trace of the velocity pressure-gradient, term is needed 
rather than a model for the pressure- transport vector. This might be easier to achieve. 


2.3.2. Fast and return splitting 
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An equation for the pressure fluctuation can be derived from the Navier-Stokes equa- 


tions 


P,kk ~ + UijUj'i u i,j u j,i 

with the Neumann boundary condition at the wall 


p'.v = v ' 


yy 


(31) 


(32) 


It is customary (see Lumley 1978, for example) to split the pressure into two parts, p' = 
p 1 -|- p 2 , one associated with the first term on the right-hand side of equation (31) and the 
other with the second and third terms. Most of the analyses used to model the pressure- 
strain terms consider homogeneous cases where the boundary conditions are not considered 
in the split. The Poisson equation and the boundary conditions are linear in p. Therefore, 
we can isolate the effects of the viscous terms at the wall by splitting the pressure into three 
parts, a “return” part, a “rapid” part, and (for the case of flows with walls) a “Stokes” 
part. ' 


i) The rapid pressure, p 1 , is defined as the solution to the following problem: 

Pjj ~ -2UiJ u 'j,i 


with the boundary conditions at the walls 


4=0 


(33 - a) 


(33 - 6) 


ii) The return pressure, p 2 , is defined as the solution to 

2 /ti } 7 \ 

Pjj = ~ 

with the boundary conditions at the walls 

4 = o 


(34 - a) 


(34 - b) 


iii) And finally the Stokes pressure, p s , is defined as the solution to 


Pji = o 


with the boundary conditions at the walls 


s _ i 
P<y ~ v *yy 


(35 - a) 


(35-6) 


This split resolves the question of whether to add the boundary conditions to the 
return part of the pressure or to the rapid part. It does not remove the effect of the wall 
on the rapid and return pressure. The pressure-strain terms are linear in p and therefore 
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the Stokes pressure-strain statistics can be added to either the rapid pressure-strain terms 
or to the return pressure-strain terms. 

The rapid part of the pressure-strain can be written analytically as follows: 

/ J J Js;., G dV (36) 

where G is the Green function with homogeneous Neumann boundary conditions at the 
walls. Note that most modelers neglect the surface integral terms that should be added 
to equation (36) if inhomogeneous Neumann conditions are used for the pressure. The 
use of homogeneous boundary conditions (eq. 33 — b) at the walls for the rapid pressure is 
consistent with equation (36) and the approximation used by the modelers. The effect of 
the wall on the pressure is contained in the form of the Green function G. 

Figures 13-16 show the splitting of the pressure-strain term into the three components. 
For the case of </> n the rapid part of the pressure-strain term is of the same order as the 
return part at y + > 80. Close to the wall the return part is larger than the rapid part. The 
<t > 22 terms show that at > 80 most of the correlation is due to the return part. Near the 
wall, the rapid part in this case has the opposite sign from the total term. The total term 
is consuming close to the wall, while the rapid part is producing. The <f> 33 terms show that 
the rapid part contributes the most to this component. Close to the wall the return part 
becomes the main contributor. The 4>\ 2 split shows that at y + > 80 the return and the 
rapid terms are of the same order. Close to the wall the return is the main contributor to 
the total term. The behavior of <f> i2 near the wall is much more complicated than that of 
<f>ij (eq. 30), which also suggests that <}>jj might be simpler to model. Except for which 
is negligible throughout the channel, for other components is significant only near the 
walls. 

3. Model testing 

In the previous sections we presented the budget data for the channel. In this section, 
we will use this data to evaluate some existing turbulence models. Our testing will be by 
direct comparisons of the terms in the budget with the model expression using the channel 
data. 

3.1. Dissipation rate models 


3.1.1. Algebraic models for e t j 

Rotta (1951a) argued that in the limit as Re »-*• 0, the dissipation rate tensor will be 
aligned with the Reynolds stress tensor and can be modeled as 


tij — C, 


U':U' 

l J 

2k 


(37) 


He also argued that in the limit Re t— ► oc, the dissipation rate tensor is isotropic. This 
idea was used by HL2 who argued that the model for the dissipation rate tensor should 
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take the following form: 


2 [ 3 u'u' 1 

€i > = r\ {1 ~ ft)6ij + 2nr f ’\ (38) 

where they inferred from the experimental data that f, is a function of the turbulence 
Reynolds number (k 2 /e) as follows: 


/. = d + ^ 7 )- (39) 

The assumed form in equation (38) implies that the anisotropy tensor of the Reynolds 
stress and the anisotropy tensor for the dissipation rate are related as follows: 

dij = fsbij (40) 

where d t j = c f j/( 2 c) - and bij = u\u'jl(‘lk) - S^/Z. 

Lumley and Newman (1977) identify a turbulence state in terms of the second (II ) 
and third (III ) invariants of the Reynolds stress anisotropy tensor {bij). They have shown 
that turbulence states can be identified on a -II vs. Ill map (anisotropy map) and that due 
to the properties of bij , turbulence states are limited inside the region bounded by the two 
axisymmetric states and the two-dimensional state (see fig. 17). It can be shown that the 
states of the dissipation anisotropy tensor (djj) are also contained in the same region as 
the bij tensor (for more detail, see Lee & Reynolds, 1985). If we neglect the d \ 3 and the d 2 3 
components of the tensor compared to the other components (for a large enough sample 
they are negligible), we can use the budget results presented in the previous section to 
compute the variation (as a function of y + ) of dij on the anisotropy map. Figure 17 shows 
the points for dij and bij on the anisotropy map for different t/ + locations in the channel. 
Figure 18 shows the points for and the curve for the right-hand side of equation (40). 
We note that the state of turbulence producing the dissipation rate tensor vary from a 
nearly isotropic state in the center of the channel to a two-dimensional state close to the 
wall. The model (eq. 40) clips the transition from the almost two-dimensional state near 
the wall to the state in the core region. Close to the center, the model is closer to the 
axisymmetric state than the data would indicate. The fact that the dissipation rate tensor 
is close to a two-dimensional state near the wall is an indication that the variation in d ij 
near the wall is due in part to wall-proximity. Near the wall, the normal component to the 
wall is suppressed and the anisotropy tensor approaches the line of the two-dimensional 
state. At around y + = 3.46 the state of the dissipation anisotropy tensor is closest to 
the one-dimensional state. The b ^ tensor will also vary from a two-dimensional state 
near the wall to the nearly isotropic state in the core region. In fact, if we assume that 
f s = 1 (i.e, that the anisotropy tensor of the dissipation rate and the anisotropy tensor of 
the Reynolds stress are equal), we find (fig. 17) better agreement between the data and 
the model. It is possible that this agreement is because KMM’s flow is at low Reynolds 
number. Comparison of the anisotropy invariant map of bij with the map of d t j shows that 
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in the core region, the dissipation anisotropy is closer to the axisymmetric state than is the 
Reynolds stress anisotropy. We point out that close to the wall, Taylor-series expansions of 
djj and b{j show that, they are equal only up to 0(t/ + ). For example, dj 2 and fc ]2 approach 
0 as y + t— ► 0, but the ratio of the two terms will yield dj 2 /6 12 i— > 2 as y t— * 0. In fact, 
Launder & Reynolds (1983) proposed a model that will have the proper limits (for the 
ratio of the component of bij and d,j ), but using the values of the constants (a and f3) 
recommended by them will yield a model that is tensorially incorrect: e.g., the trace of tij 
and the trace of their model are not equal. 

Figures 19-22 show the four components of £j j compared to the components of the 
model, — e u'ju'j/k. The off-diagonal component shows the largest difference between 
the model and the data. The diagonal components show better agreements but would 
require a different function, /,, for the different components to obtain an improvement in 
the agreement. 

3.1.2. Transport models for e. 

Almost any type of one point closure models would require a time-scale or a length- 
scale model. Often, the dissipation rate of the turbulence kinetic energy is used to obtain 
these scales. In addition to the equations for the transport of the Reynolds stresses, an 
equation for the transport of the trace of e; j which is twice the dissipation rate of the 
turbulence kinetic energy, e (see eq. 23) is used. The terms in the equation of c have been 
modeled by a number of workers (see for example Davydov, 1961; DH; HL1; Lumley & 
Khajeh-Nouri, 1974; HL2). Most of the current models for the e-equation can be written 
as the sum of a production term, a dissipation term, a turbulent transport, and a viscous 
diffusion term 

- C«J + C, (^ c,,) + €jj (41) 

where P & = Pu/2 is the production rate of the turbulence kinetic energy, C t j, C f2 , and 
C ( are constants or functions of the turbulence Reynolds number. The disagreement is 
in the correspondence of the modeled expression with the exact equation. Lumley and 
Khajeh-Nouri (1974); and later HL2, associated the right-hand side of equation (41) with 
the difference between P* and T (see eq. 23) and did not identify a model for each term 
in the balance equation. Davydov (1961) and HL1 combined the production terms P} and 
P 2 and modeled them as follows: 

p; + p; = -C„ = c„ l Pi, (42) 

HL2 recommend C (1 = 1.275. Figure 23 shows the above model compared to the data; 
the agreement is good away from the wall. The model predicts the production very well 
away from the wall. In the near wall region, however, the peak in the production of e is 
underpredicted and a modification to the model is needed in this region. 

HL1 combined the production term P* with T and modeled the combination as 

-P?+r = C (2 j ( 43 ) 
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The right-hand side of equation (43) approaches oo as y *-> 0. HL2 modified the model 
for the dissipation rate of e using a modified dissipation rate of k, e = e - 2 ((fc 1 / 2 ) y +) , 
which insures that the ratio l/k is bounded as y + i— *■ 0. They also argued that turbulence 
data indicate that the dissipation rate of e is a function of the Reynolds number. This 
functional was accounted for by introducing a damping function. If the model of HL2 
represents a model for the left-hand side of equation (43), we have 

-P < 4 + T = C (2 /,| (44) 

where C ( 2 = 1.8, and f ( = 1. — exp [— (fc 2 /6e) 2 ] . Figure 24 shows the above model 
compared to the data. In the core region, the model and the data show good agreement. 
Close to the wall, the model underpredicts the data. 

The turbulent transport of e is modeled by HL1 as 



where C t = 0.15. Figure 25 shows the model compared to the data. This term is small 
compared to the other terms in the budget equation and the disagreement between the 
model and the data is small compared to the errors in the other terms. The production 
rate, P 3 , and the pressure diffusion rate, II t , were neglected by HLl, and the present data 
also show that these terms are small. 

3.2. Pressure-strain models 


Most models used for the velocity pressure-gradient expression are based on splitting 
the expression into a pressure-strain term and a pressure diffusion term. The pressure 
diffusion term is either added to the turbulent transport term or neglected. Several models 
for the pressure-strain term exist that use different approximations and arguments to 
provide closures for the term. Most of these closures are based on homogeneous flow 
arguments. In this section, we will test the closure of LRR for the pressure-strain term 
that was developed for wall-bounded flows. LRR split the pressure-strain term into a 
return term, a rapid term, and a wall term. For the return terms, c> 2 ; . they recommend 
the use of the model proposed by Rotta (1951a), 


= ( u Wj - (46) 

where (7j is a model constant. They modeled the rapid term by assuming that, the mean 
velocity gradient is slowly varying and write 

t'i = (47) 


They then assumed that a™ is linear in the Reynolds stresses. Substitution of the linear 
approximation for a ™ into the expression for (f>\- yields 



(C 2 + 8) 

11 


Pij 3 ^^ 

(SOO, - 2) _ (80,-2) 

55 1 *’ J 1 11 

2 

■^ij ~ g Pk^ij 



(48) 
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where, Aij = —{u\u' k Uk,j + u'ju' k Uk,i). The above model also has one adjustable constant, 
C 2 . The value of the constants C\ =1.5 (for the return term) and C 2 — 0.4 were chosen 
by LRR by matching the homogeneous shear experiment of Champagne, Harris, & Corrsin 
(1970). 

By examining the case of wall bounded flows, LRR argued that a third term, 
is needed to account for near-wall effects, corresponding to the reflected wall influence of 
4>}j + 4>\ j. They argued that the model for the wall effects should take the same form as 
and using near-wall data, the wall effect on the pressure-strain was modeled as 



k 3 * 

‘ 1 1 

<t>?j = 1 0.125- U'u' - -kSij) +0.015 (Pij - Aij) 

It 

[y+ + (2Re r -y+)\ 


Figures 26-28 show comparisons of LRR’s model for the diagonal terms to the data. 
The comparison for the off-diagonal term is similar, the agreement is acceptable away 
from the wall, but is poor close to the wall. We have also shown the distribution of the 
individual terms. Rotta’s return model does not vanish at the wall; while, as can be seen 
from the Taylor series expansions (eq. 26), the diagonal terms should vanish at the wall. 
This is an indication that LRR’s model will behave poorly close to the wall because of 
the return model. It is clear that Rotta’s model should be modified to include the correct 
behavior near the wall. In fact, the model without Rotta’s return model shows the proper 
trends, and it will be better than the full model. 

3.3. Turbulence transport models 

The nonlinearity of the equations of motion introduces higher order moments when 
equations for the moments are derived. For the Reynolds-stress equations, the triple- 
correlation terms need either to be closed or to have an equation derived for the m. The 
need for equations describing the evolution of the turbulent transport term ( Tijk = t^tt'V fc ) 
was suggested by Chou (1945a) and Davydov (1959). But using closures at the triple- 
correlation level will add ten more equations to the system of equations to be solved. How- 
ever, the transport equations for Tijk are used to derive closures for the triple-correlation 
terms. The construction of the model for the transport terms starts (see HL1; and Lumley, 
1978) with the governing equation for the transport terms, 


D 

Dt 


Tijk + {TijiUk.i + T k iiUj,i + TjkiUij} - {«'«' («),/ + + «><(“X).'} 


+ (u\u , i u'jU'k),l = ~{ U 'j U 'kP[i + U 'i U 'jP[k + U 'k U 'iP[j } + K«// + U i U 'j U 'k,U + U 'k U i U 'j,ll} 

(50) 


Davydov neglected the second term on the right-hand side of equation (50). He also 
argued that in an analogy to the pressure-strain model, the remaining term on the right- 
hand side should be modeled as —C' D e/k Tijk- To close the quadruple-correlation term, 
Millionshchikov’s zero- fourth- cumulant hypothesis (Monin & Yaglom 1975, p. 241) is often 
invoked 


= u i u j u 'k u 'i + u i u 'k u 'j u \ + u i u 'i u 'j u 'k 


(51) 
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Using the above closure in the triple-correlation equations, we have 


D 


jjfTijk + {T t jlUk,l + TkilUj'l + TjklUij} + {u' i u'i(u , jU , k ) t i + u'jU'i(u'i t U , i ) t l + u' k u , l (u' i u'j) j} 




(52) 


If we use the assumption of HLl (that the transport terms are in equilibrium) and 
drop the term, the model for the triple correlation will close as follows: 

+ PwJTyi - {jCojiii + U,j}T, kl - {icl.jiu + Uj,,}T til 

= {Ku\{u , j u' k ) il + u' k u' l (u' i u' j )' l + ttjuj(tt^) f i} 

In addition to the equilibrium assumption, HLl assumed that the production terms are 
negligible and wrote the model for the triple correlation as follows 




(54) 


For the case of no mean velocity gradients, the two models (eqs. (53) and (54)) are the 
same (with C' D = 1/(7,). If we use the expression given by equation (53) to model the 
transport terms, we have for the channel case 


t 2U = c s - 
£ 


k r 


2u\u' 2 (u'jii' 2 ) )2 + u' 2 u' 2 (uX),2 + 2C/ 1 , 2 7 , 212 


-T 222 = C t -u' 2 u' 2 (u^^), 2 


fV 

—T233 — 


k 

e L 


2 U 2 u’ 3 (u 2 « 3 ), 2 + u' 2 u' 2 (u' 3 u ' 3 )'2 


-T212 = c.~ 


k r 


£ L 


2 u 2 u 2 (uj u 2 ) t 2 + u\n' 2 {u' 2 u 2 )'2 -I- u J)2 T 222 


The above expressions show that the production of the triple correlation in the fully de- 
veloped channel will affect the T 2 n and T 212 components only; the largest effect is on the 
r 2 ,i component. Figure 29 shows the models for r 2 n given by equations (53) and (54) 
(with C s = l/C' D = .11) compared to the data. These results indicate that both models do 
not agree well with the data and that including the production term does not improve the 
model. The extra effort involved in inverting the coupled system given by equation (53) is 
not justified. 


A simpler form for the trans port terms was derived by DH using the recipe that 
turbulence transport of a quantity, should be modeled as cx u' k u\ <j> t . Following this 

recipe, they modeled the transport term as 



Although this model is tensorially incorrect, as pointed out by HL1, it is often used 
because of its simplicity. LRR recommend C a = .11 and C' a = .25 (based on numerical 
experimentation) for the values of the constants in the models. Figures 30-33 show the 
comparisons of the models given by equations (54) and (55) with the data. We can see 
that the simpler (and tensorially incorrect) model of DH performs as well as HLl’s model 
and, for the case of the T-ju component, better than the model of HL1. Considering the 
simplicity of both models, the trends predicted by the model are acceptable except for T 033 
where the sign of the slope close to the wall is not predicted correctly. 

We note that both HLl’s and DH’s models do not have the proper behavior near a 
wall; for example u' 2 u' 2 u 2 = 0(y +6 ) as 1 — ► 0. The different terms in the expressions of 

the model will asymptote as y + 1 — > 0 as follows: 

k=0(y+ 2 ) 

u 2 ii^ =0(y +3 ) for / / 2 
u' 2 u ' 2 = 0 (y +4 ) 

Therefore, 

~{u' 2 u 2 u' 2 u 2 ) t2 - 0 (y +9 ) 

The model for u 2 u 2 u ' 2 will not have the proper behavior as y + 1 — ► 0. However, the asymp- 
totic behavior only holds extremely close to the wall ( y + < 5) and should not be the sole 
determining factor in evaluating a model. 

4. Summary and discussions 

We have used the channel data of Kim, Moin & Moser (1987) to compute the terms 
in the transport equations for the Reynolds stresses and to compute the terms in the 
transport equation for the dissipation rate of turbulence kinetic energy. It is important 
to recognize that the simulation data is at low-Reynolds number and that the present 
conclusions maybe valid only for low-Reynolds number flows. In particular, the results 
near the wall related to dissipation rates are sensitive to Reynolds number dependence. 
The budgets reveal that several terms that were negligible away from the wall become 
important close to the wall. The budget for the turbulence kinetic energy reveals that, 
contrary to the commonly held belief (see for example Townsend, 1976), the pressure-gain 
term near the wall is not of the same order as the production term away from the wall. 
In fact, the pressure term remains small relative to the dissipation rate and the viscous- 
diffusion rate. Away from the wall, the budget for the dissipation rate of the turbulence 
kinetic energy reveals that the turbulent production term and the dissipation rate are the 
dominant terms as estimated by Tennekes & Lumley (1972). Close to the wall, all the 
production terms become important. In considering the velocity pressure-gradient term, 
we have shown that different splits of the velocity pressure-gradient term will give different 
behaviors near the wall and that a judicious choice may be necessary. 
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For wall-bounded flows, we have shown that the inhomogeneous boundary condition 
on the pressure introduces a third term in the split of the pressure and have recommended 
that the pressure be split into a rapid term, a return term, and a Stokes term. We find 
that the rapid and the return terms in the channel are of the same order, and we cannot 
neglect one with respect to the other. The 22 component of the pressure-strain term shows 
that the rapid part in fact has the opposite sign as the total term. Away from the wall 
the rapid part is consuming, while the total term is producing. Close to the wall the total 
term is consuming (the splatting effect), while the rapid part is producing. 

Comparison of closure models with the data reveals that the pressure-strain term 
needs immediate attention and the model of Launder, Reece & Rodi (1975) has difficulty. 
As a first approximation, the anisotropy tensor for the dissipation rate of the Reynolds 
stresses may be modeled in terms of the anisotropy tensor of the Reynolds stresses. The 
budget for the dissipation rate of the turbulence kinetic energy is modeled well away from 
the wall; close to the wall, improvements are needed. Finally, the transport term also can 
be improved upon. Overall, the closure models are better than expected for the budget of 
the dissipation rate of the turbulence kinetic energy and are generally inadequate for the 
pressure-strain correlations. 
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Figure 1 . Terms in the budget of u\u\ in wall coordinates. P n = Production; 
Tn = Turbulent transport; £>n = Viscous diffusion; en = Dissipation rate; Hu = Velocity 
pressure-gradient term. 



Figure 2 . Terms in the budget of in wall coordinates. T22 = Turbulent 

transport; D22 = Viscous diffusion; e 2 2 = Dissipation rate; II22 = Velocity pressure- 
gradient term. 
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Figure 3. Terms in the budget of u' 3 u 3 in wall coordinates. T 33 = Turbulent 
transport; D 33 = Viscous diffusion; e 33 = Dissipation rate; II 33 = Velocity pressure- 
gradient term. 



Figure 4. Terms in the budget of in wall coordinates. P J2 = Production; 
Tj 2 = Turbulent transport; D 12 = Viscous diffusion; e 12 = Dissipation rate; IIj 2 = Velocity 
pressure-gradient rate. 
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Figure 5. Terms in the budget of the turbulence kinetic energy, k , in wall co- 
ordinates. Pk = Production; 7* = Turbulent transport; Dk = Viscous diffusion; tk = 
Dissipation rate; II & = Velocity pressure-gradient term. 



Figure 6. Terms in the budget of the dissipation rate of the turbulence kinetic 
energy, e, in wall coordinates. P} = Production by mean velocity gradient; P 2 = Mixed 
production; P { 3 = Gradient production; P* = Turbulent production; T ( = Turbulent 
pransport; D t = Viscous diffusion; T = Dissipation rate; !!« = Pressure transport. 
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Figure 7. Split of the viscous term into a dissipation rate term, T, and a viscous 
diffusion term, D ( . 



Figure 8. Split of the velocity pr essure- g radient term into a pressure transport 
term, (pu ' 2 ) t „+, and a pressure-strain term, pu' + . 
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Figure 9. Split of the velocity pressure-gradie nt ter m, Du, into a pressure transport 
term, u\u\/2k IIu, and a redistributive term, IIu — u' i u\/2k II». 



Figure 10. Split of the velocity pressure-gradient term, II 22 , into a pressure 

transport term, u' 2 u' 2 /2k II jj, and a redistributive term, II 22 — u 2 u 2 /2k IIu . 
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Figure 11. Split of the velocity pressure-gradient term, II33, into a pressure 

transport term, u' 3 u' 3 / 2 k II//, and a redistributive term, II33 — u 3 u 3 / 2 k II//. 



Figure 12 . 


Split of the velocity pressure- gradient term, II j 2 , into a pressure 


transport term, u' 1 u' 2 / 2 k II//, and a redistributive term, n 12 — u\ u' 2 / 2 k II//. 
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Figure 13. Split of pressure-strain term, <j>n, into a rapid term, <j>\ 
and a Stokes term, <f> f,. 



Figure 14. Split of pressure-strain term, <f > 221 into a rapid term, <f>. 
<j> 2 2 , and a Stokes term, </>f 2 - 




Figure 15. Split of pressure-strain term, <f> 33 , into a rapid term, <f> 33 , a return term 
<f>l 3 , and a Stokes term, <f> 33 . 
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Figure 16. Split of pressure-strain term, </>i 2 , into a rapid term, <f>] 2 , a return term 
4>] 2 , and a Stokes term, <f) 32 . 




Figure 17. Anisotropy invariant map. o o o o dij at various y + in the channel; 
- bij at various y + in the channel. 
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Figure 18. Anisotropy invariant map. o o o o dij at various y + in the channel; 
- model, equation (40). 
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Figure 19. Distribution of en across the channel, ooooeu term computed from 
the channel data; model, e/k u\u\. 
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Figure 20. Distribution of e 22 across the channel, o o o o c 22 term computed from 
the channel data; model, e/k uLul. 




Figure 21. Distribution of €33 across the channel. 0000 633 term computed from 
the channel data; model, e/k u' 3 u' 3 . 
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Figure 22. Distribution of e 12 ac ross the channel, o o o o e 12 term computed from 
the channel data; model, e/k u\u' 2 . 
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Figure 23. Distribution of the production term, P} 4- P? , in the budget of e across 
the channel; o o o o term computed from the channel data; model, equation (42). 



Figure 24. Distribution of the dissipation term, P '* — T, in the budget of e across 
the channel, o o o o term computed from the channel data; model, equation (44). 
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Figure 25. Distribution of the turbulent transport term, T ( , in the budget of e 
across the channel, o o o o term computed from the channel data; model, equa- 

tion (45). 



Figure 26. Pressure-strain term, pu[ j, in the budget equation for ttjUj across the 
channel, o o o o term computed from the channel data; 

model, (eq. (46)+eq. (48)+eq. (49)); model, equation (46); model, 

(eq. (48)+eq. (49)). 
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Figure 27. Pressure strain term, pu' 2 2 , in the budget equation for u' 2 u' 2 across the 
channel, o o o o term computed from the channel data; 

model, (eq. (46)+eq. (48)+eq. (49)); model, equation (46). model, 

(eq. (48)+eq. (49)). 



Figure 28. Pressure strain term, pu' 3 3 , in the budget equation for u' 3 u 3 across the 
channel, o o o o term computed from the channel data; 

model, (eq. (46)+eq. (48)+eq. (49)); model, equation (46). model, 

(eq. (48)+eq. (49)). 
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Figure 29. Triple correlation term -u' 2 ttj«i across the channel. 

o o o o term computed from the channel data; model, equation (53); model, 

equation (54). 



Figure 30. Triple correlation term —u' 2 u\u\ across the channel. 

o o o o term computed from the channel data; model, equation (54); model, 

equation (55). 
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Figure 31. Triple correlation term — u' 2 u' 2 u' 2 across the channel. 

o o o o term computed from the channel data; model, equation (54); model, 

equation (55). 



Figure 32. Triple correlation term —u' 2 u' 3 u' 3 across the channel. 

oooo term computed from the channel data; model, equation (54); model, 

equation (55). 
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Figure 33. Triple correlation term —u^u^ u' 2 across the channel 

oooo term computed from the channel data; model, equation (54); m 

equation (55). 
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